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Abstract: The D2 variable region of 28S ribosomal RNA was sequenced from ethanol preserved specimens or obtained from the literature to pro- 
vide the first molecular phylogenetic reconstruction of this subfamily. Evolutionary relationships were investigated by comparing the results using 
three different analysis methods < distance-based neighbor-joining, NJ; maximum parsimony, MP; and maximum likelihood, ML) and three differ- 
ent outgroups. The monophyly of the Rogadinae s. s. is well supported by all trees generated from molecular data. The morphologically defined 
tribes and the relationships within the subfamily are partially resolved. We found no evidence for a monophyletic Rogadini lineage, neither for the 
morphological subtribe division within the tribe, but we did find some evidence for three other tribes. Among the genera of Rogadini our results 
show strong support for the paraphyly of this group. It can be divided into two clades, one consisting of Aleiodes + Arcaleiodes, the other of 
Canalirogas + Conspinaria + Spinaria + Rogas. The clade Aleiodes + Arcaleiodes was constantly resolved using different outgroups with all 
Aleiodes species forming a monophyletic group. Two distinct clades are found within Aleiodes» strongly corresponding to the morphologically de- 
fined subgeneric divisions of the genus. For the clade Canaliregas + Conspinaria + Spinaria + Rogas, we found only limited support for this clade 
to be a monophyletic group. We did not find any support for the genus Spinaria to receive subtribe rank. The tribe Clinocentrini is well supported 
in our analyses, and probably occupy a basal position in the phylogeny of the subfamily. We also found two well-supported clades containing the 
genera Yelicones and Stiropius, corresponding to the tribes Yeliconini and Stiropiini, respectively, which were well defined by morphological data, 
but the placement of these two clades within the phylogeny of the subfamily is poorly resolved in this study. 
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considered to be monophyletic (van Achterberg» 1991) and 


1 Introduction 


The Rogadinae sensu lato is a very large subfamily, 
once including a number of cyclostome tribes that differ widely 
in biology and share few synapomorphies with one-another 
CShaw. 1983: van Achterbergy 1984: Shaw and Hud- 
dleston> 19912. It therefore has recently been considered a 
heterogenous group of polyphyletic C Quicke and Achterberg» 
1990: van Achterberg; 1991), or paraphyletic C Shaw and 
Huddleston, 1991), nature. The concept of Rogadinae sensu 
stricto L = Rogadini of most authors | is widely accepted at pr- 
esent; and defined biologically as all those braconid species 
that mummify lepidopteran host larvae (Whitfield, 1992: van 
Achterberg; 1993: Shaw. 1997: Chen and He; 1997). It is 


furthermore was supported by a phylogenetic analysis of rela- 
tionships in the cyclostome lineages of the Braconidae by 
Whitfield (19922. However: the phylogenetic relationships 
among the genera and the evolution of the group in general re- 
main poorly understood. The most current classification 


scheme for Rogadinae s. s. is that of van Achterberg 


(1993), who divided the subfamily into four tribes: i. e.> 
1991. Rogadini Foerster, 
1862, Stiropiini van Achterberg» 1993, and Yeliconini van 


Clinocentrini van Achterberg» 


Achterberg, 1991. Rogadini is the largest of the four tribes 
and includes the majority of known genera and species: it is 
subdivided into two subtribes: Spinariina and Rogadina. 
Because of its high diversity, estimating a phylogeny of 
this subfamily is a complex task. To date no molecular phylo- 
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genetic study has been concentrated at the generic level in 
this group. 

In this paper we use the 285 gene to resolve the phylo- 
genetic relationships within this group of parasitic wasps. The 
285 gene appears to be a promising molecular maker that has 
been useful for estimating relationships corresponding to ge- 
neric and subfamily levels of divergence within the family Bra- 
conidae (Belshaw and Quicke; 1997a; b: Dowton and Aus- 
tin; 1998; Mardulyn and Whitfield, 1999). 


2 Materials and Methods 


2.1 Taxa examined 

We examined 28 species belonging to 17 genera. The 
species and their provenance are listed in Table 1. Identifica- 
tions of the taxa for gene sequencing were made by the first 
author. The tribe and subtribe arrangement followed in this 
paper is that of van Achterberg (1993) and Chen and He 
(1997) based on morphological and biological characters. 
Our analysis therefore includes at least one representative of 
all the four tribes and two subtribes of Rogadini recognized by 
van Achterberg (1993) and Chen and He (1997). 

DNA sequences in this paper have been deposited in 
EMBL/GenBank/DDB] database under accession numbers 
listed in Table 1. 

2.2 DNA amplification and sequencing 

Sequences were obtained either from the literature C see 
Table 1 for sources) or generated from the taxa that appear in 
Table 1. Genomic DNA was extracted from insects stored in 
100% ethanol. Legs were removed from larger wasps È > 5 
mm long) and used for exactions whereas the abdomen and 
thorax were used for smaller wasps € < 5 mm long). Ethanol 
was then removed by washing 3 times (15 min each) in 10 
mmol/L Tris-HCl C pH 8.0) containing 100 mmol/L NaCl and 
1 mmol/L MgCl, . Tissue was then ground in 400 pL of 10 
mmol/L Tris-HCl (pH 8.0); 10 mmol/L EDTA; 1% SDS; 
and then 100 pg Cor 6 uL) proteinase K and incubated at 
37°C for 1 h. The homogenate was extracted with phenol: 
chloroform: isoamyl alcohol (25:24:12. DNA was redissolved 
in 100 pL of sterile water and stored at — 20°C. 

A portion of the 285 rDNA gene was then amplified. 
PCR reactions were carried out in 50 uL volumes containing 
100 ng extracted DNA, 5 pL 10 x Taq buffer; 4 pL 25 
mmol/L MgCl,» 4 uL 2.5 mmol/L dNTPs, 1 uL 20 pmol/L 


primers, 0.5 uL Taq DNA polymerase (5 U/L). The prim- 
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ers used for all amplification followed Belshaw and Quicke 
(1997): 5’-AGA GAG AGT TCA AGA GTA CGT G-3’ Cor- 
ward primer) and 5'-TTG TGC CGT GTT TCA AGA CGG G- 
3’ Creverse primer). PCR conditions were: 30 cycles of 94°C 
denaturation (30 s); 54°C annealing (1 min); and 72°C ex- 
tension (1 min) with an initial denaturation of 3 min at 94°C 
and a final extension for 8 min. PCR products were purified 
with the QIAquick PCR purification kit and then sequenced 
directly in both directions using the same primers and TaqFS 
Dye-labeled terminator on an ABI 377 automated sequencer. 
2.3 Sequence alignment 

New sequences were deposited in GenBank under acces- 
sion numbers from AY167645 to AY167659. Sequences were 
aligned based on distance approach using CLUSTAL X version 
1.81 CThompson ef al.» 1997) with different gap opening 
and gap extension values. Because some sequences used in 
this paper from the literature contain both regions D2 and D3 
or other regions (D2 — D10 of Apis mellifera)» manual align- 
ment was used to remove the D3 or other regions. Lengths of 
the resulting alignments were between 499 to 560 bp. 
2.4 Phylogenetic analysis 

Following alignment, three different methods of phyloge- 
netic analysis were performed using PAUP* 4.0 Cbeta 10 ver- 
sion Swofford, 20012. First; we used maximum parsimony 
(MP). To find the most parsimonious tree (s); heuristic 
parsiomony searches (Hillis et al.» 1996) were performed 
using 100 replicates of random addition sequences and the 
TBR option for branch swapping. Each base was treated as an 
unordered character with equal weights, with gaps treated as 
missing data. Where more than one most parsimonious tree 
was found, a strict consensus tree was calculated. Down- 
weighting transitions or treating gaps as a fifth base did not 
markedly affect the results obtained. Statistical support for 
each node was evaluated by bootstrap analysis ( Felsenstein» 
1985) with 1 000 replicates. Second, we used a distance- 
based method following the neighbor-joining algorithm (NJ) 
(Saitou and Nei» 1987) for obtaining a minimum-evolution 
tree. Bootstrapping evaluation of each node was performed as 
above. Thirds we employed maximum likelihood (ML). ML 
trees were generated under the HKY85 model using base fre- 
quencies estimated by PAUP: and the default number of sub- 
stitution types (2, HKY85 variant) and transition/transversion 
ratio (2). Heuristic searches were used with 10 replicates of 
random addition sequences and TBR branch swapping. Boot- 


strap analysis was performed with 100 replicates. 
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Taxa name 

Rogadinae 

Clinocentrini 
Clinocentrus politus Chen et He 
Clinocentrus sp. 

Ydiconini 
Yelicones wui Chen et He 
Yelicones sp. 

Stiropiini 
Stiropius sp. 

Rogadini 

Spinarina 
Spinaria fuscipennis Brullé 
Rogadina 

Aleiodes aestuosus C Reinhard) 
Aleiodes convexus van Áchterberg 
Aleiodes cariniventris (Enderlein) 
Aleiodes miroculatus (Watanabe) 
Aleiodes altemator (Nees) 
Aleiodes coxalis (Spinola) 
Aleiodes esenbecki Hartig 
Aleiodes gracilipes (Telenga) 
Aleiodes mythimnae He e Chen 
Aleiodes similes ( Curtis) 
Arcaleiodes aglaurus (Chen et He) 
Conspinaria sp. 
Rogas sp. 
Canalirogas sp. 

Braconinae 
Amyosoma. chinensis (Szépligeti) 
Vipio sp. 

Opiinae 
Opius caricivorae Fischer 

Histerominae 
Histeromerus mystacinus 

Rhyssalinae 
Rhyssalus sp. 

Ichneumonidae 
Agrypon japonicum Uchida 
Itoplectis naranyae (Ashmead) 

Apidae 
Apis mellifera L. 


based on the D2 variable region of 28S ribosomal RNA 


Table 1 Taxa examined 


Abbreviation in Figs. 1 — 3 (collected locality) 


Clinocenpo (Zhejiang, China) 
Clinocensp 


Yeliconesw (Guizhou, China) 


Yeliconess 


Stiropiuss 


Spinariafu 


Aleioaestu Guizhou; China) 
Aleioconve (Guizhous China) 
Aleiocarin (Guizhou, China) 
Aleiomiroc (Zhejiang: China) 
Aleioalter 

Aleiocoxal 

Aleioesenb (Shandong, China) 
Aleiograci (Guizhou, China) 
Aleiomythi (Zhejiang: China) 
Aleiosimil 

Arcaleioag (Zhejiang. China) 
Conspinari 

Rogassp 

Canaliroga 


Amyosomach (Zhejiang: China) 


Vipiosp (Inner Mongolia. China) 


Opiussp (Zhejiang, China) 


Histeromer 


Rhyssaluss 


Agryponjap (Zhejiang: China? 
Itoplectis (Zhejiang; China? 


Apismellif 
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Accession number 


AY167658 
[ AJ302890 |* 


AY167659 
[ AJ302905 ]* 


[ AJ302904 }* 


[797064 P 


AY167655 
AY167654 
AY167650 
AY167653 
[ AJ302884 )* 
[283611 I 
AY167651 
AY167652 
AY167656 
[ AJ302885 )* 
AY167657 
[ AJ302898 }* 
[AJ296058 I 
[ AJ302889 }* 


AY167647 
AY167648 


AY167649 


[ 283601 |! 


[ 283603 |! 


AY167645 
AY167646 


[C AJ302936 * 


1. Sequences from Belshaw and Quicke, 1997; 2. Sequences from Belshaw e al., 1998; 3. Sequences from Belshaw & al., 2001; 4. Sequences from 


Belshaw and Quicke: 2002. 
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2.5 QOutgroup selection 

We used 3 different outgroups to analyze tree root- 
ing. The most distant groups were two icheumonid species 
CAgrypon japonicum, and ltoplectis naranyae ) and one 
aculeate species C Apis mellifera). The Ichneumonidae is 
widely accepted as the monophyletic sister-taxon to Bra- 
conidae both by morphological data CSharkey and Wahl. 
1992; Gauld and Shaw. 1995) and by molecular evi- 
dence CDowton and Austin. 19942. and the aculeate of 
Hymenoptera is considered as the probable sister-group of 
the Ichneumonoidea C Rasnitsyn. 1988; Sharkey and 
Wahl. 1992; Gauld and Shaw, 19952. The second group 
was the Rhyssalinae and Histeromerinae: two subfamilies 
of Braconidae that are considered to occupy a basal posi- 
tion within the family © Quicke and van Achterberg» 
1990; Whitfield. 19922. The third was the Braconinae 
C Amyosoma chinensis» and Vipio sp.) and Opiinae 
COpius caricivorae Fischer); both are relatively close to 


Rogadinae within the monophyletic cyclostome clade. 


3 Results and Discussion 


We tested alignments using the Clustal X program with 
different gap opening and gap extension values: and resulted 
in different lengths of aligned sequences. The result was 
identical to that of Morrison and Ellis (1997). They conclud- 
ed in their paper that the multiple alignments from different 
procedures varied greatly in length and the alignments pro- 
duced from the Clustal W program with different gap weights 
were at least as different from each other as those produced by 
the different alignment algorithms C Morrison and Ellis, 
1997). Because the default parameters in version 1.81 Cgap 
opening 15, gap extension 6.66) were optimized using the 
balibase multiple alignment benchmark and gave the best 
overall results in the 142 alignment tests in balibase (J. 
Thompson, pers. comm.» we used the alignments with de- 
fault parameters for analysis in the present paper. 

The D2 region of 285 rDNA amplified from different taxa 
ranged from 449 to 482 bp in length: the aligned sequences 
were different in length when different outgroups were used. 
For outgroup one; the aligned sequence was 540 bp due to 
introduction of gaps among sequences. Of 540 total characters 
analyzed. 224 characters (41.5%) are constant» and 206 
characters are parsimony informative (38.2%). Regarding 


the base composition: the overall CG content is 44.6946 ; 
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ranging from 38.16 % to 49.59 96; similar to the results of 
Belshaw et al. (1998): 43% in Braconidae. The base com- 
position of this region is known to vary widely among other or- 
ganisms> with CG contents of between 26% and 85% record- 
ed C Hancock and Dover, 1988). 

The trees resulting from PAUP* analyses are presented 
in Figs. 1-3. We also show on these figures the bootstrap 
values obtained from the identical analysis. 

From all trees generated from molecular data using dif- 
ferent analysis strategies and different taxa as outgroups for 
comparsion: the monophyletic nature of the Rogadinae s. s. 
is well supported; suggesting that the present defined subfam- 
ily Rogadinae is a natural group (Figs. 1 - 3). 

The morphological-defined tribes and relationship within 
the subfamily are partially resolved. We found no evidence 
for a monophyletic Rogadini lineage: nor for the morphologi- 
cal subtribe division within the tribe: but we did find some 
evidence for the three other tribes. 

Among the genera of Rogadini our results show strong 
support for a paraphyletic nature of this group. It can be 
roughly divided into two clades: one consisting of Aleiodes + 
Arcaleiodes: the other of Canalirogas + Conspinaria + Spina- 
ria + Rogas» both of which may be given tribal rank when 
more evidence is available to support this branch in the fu- 
ture. 

The clade Aleiodes + Arcaleiodes was constantly resolved 
using different outgroups (Figs. 1-3), suggesting its mono- 
phyletic nature. All Aleiodes species also form a monophyletic 
group with 70 — 96 bootstrap support. Within Aleiodes» two 
clades are obvious» one containing six species and the other 4 
species, with 94 — 100 and 71 — 94 bootstrap value support; 
respectively. This result strongly recovers the morphological- 
ly-defined subgeneric divisions of the genus (van Achterberg» 
1991; Chen and He: 1997: Belokobylskij and Tobias: 
2000). We also found evidence from the molecular data for 
Arcaleiodes Chen and He: 1997 to be a monophyletic and 
distinct genus (Figs. 1 — 32. 

For the clade Canalirogas + Conspinaria + Spinaria + 
Rogas» we found obvious support for this clade to be a mono- 
phyletic group (Figs. 1A» 1B, 2A, 2B; 3A; and 3B, with 
64 — 90 bootstrap support value), but we did not find any 
support for the genus Spinaria to be accorded subtribe rank . 
The latterwas first given tribal rank (van Achterberg» 1988: 
Quicke and van Achterberg» 1990) due to its remarkably dif- 


ferent morphological characters ( metasomal tergites form a 


253] CHEN Xue-Xin et al.: A molecular phylogeny of the subfamily Rogadinae ( Hymenoptera: Braconidae) 
based on the D2 variable region of 28S ribosomal RNA 213 


carapace); then sunk to subtribal rank (van Achterberg» (Fig. 1B, 1C) when Ichneumonidae and Apidae» and Rhys- 


1991). The relationships between these genera are still un- salinae and Histeomerinae were used as outgroups: but when 
certain. We found their phylogenetic relationship as Rogas + the Braconidae and Opiinae were used as outgroups: the rela- 
Spinaria + © Conspinari + Canalirogas) (Fig. 1A» 2B) or tionships became Canalirogas + Rogas + ( Conspinaria + 


Rogas + (© Spinaria + C Canalirogas + Conspinari)) (Fig. Spinaria) (Fig. 3A) or Canalirogas + © Rogas + © Con- 
2A) or Spinaria + ( Rogas + ( Conspinaria + Canalirogas )) spinaria + Spinaria)) (Fig. 3B). 
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Fig. 1 Phylogenetic analysis using NJ (A), MP (B) and ML CC) methods with Ichneumonidae and Apidae as outgroups 
Numbers above branches are bootstrap values ( 96) 
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Fig. 2 Phylogenetic analysis using NJ (A), MP CB) and 


ML (C) methods with Rhyssalinae and Histeromerinae as outgroups 


Numbers above branches are bootstrap values ( 96) 
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Fig. 3 Phylogenetic analysis using NJ (A), MP (B) and ML CC) methods with Braconinae and Opiinae as outgroups 
Numbers above branches are bootstrap values ( 962 
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The phylogenetic relationships of both the clade Aleiodes 
+ Arcaleiodes and the clade Canalirogas + Conspinaria + 
Spinaria + Rogas with other clades (tribes) within this sub- 
family are largely unresolved, although the clade Aleiodes + 
Arcaleiodes always occupied a terminal position in the phylog- 
eny of the subfamily in most of our analyses (Figs. 1 — 3). 

The tribe Clinocentrini is well supported by our analy- 
ses. It occupied a basal position in the phylogeny of the sub- 
family when Ichneumonidae and Apidae was used as outgroup 
(Fig. 1: A: B. C), corresponding to the placement of 
Whitfield (1992). However: the phylogenetic relationships of 
this clade with other groups within the subfamily are poorly 
resolved in our other analysis; even becoming the sister- 
group to the clade Aleiodes + Arcaleiodes (Fig. 3A), a re- 
sult never found from morphological characters . 

We found a well-supported clade containing the genus 
Yelicones», corresponding to the tribe Yeliconini; which is 
welldefined by morphological data Ca group with the tarsi 
distinctly shortened). Although the placement of this clade 
within the phylogeny of the subfamily is poorly resolved: we 
did find that it was a sister-group to the clade Rogas + 
CSpinaria + © Canalirogas + Conspinaria 2) (Fig. 2A) or 
sistergroup to Stiropius + © Canalirogas + Conspinaria + 
Spinaria + Rogas) (Figs. 1B» 1C). 

The tribe Stiropiini» a morphologically and biologically 

well-defined group: is also well resolved in our analysis. But 
its placement within the subfamily is largely unclear: as is its 
situation with the tribe Yeliconini. It is interesting to point out 
that in one analysis it was the sister-group of the clade 
Canalirogas + Conspinaria + Spinaria + Rogas (Figs. 1B; 
1C); suggesting the similar result presented by Whitfield 
(1992) that Stiropius-group is a sister group of Rogadini s. 
s. But this group has been hypothesized to occupy a basal 
position within the phylogeny of the subfamily (van Achter- 
berg; 1991). 
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fi: RA FAIRY 28S rRNA D2 FERIA Fiber P3 REIP E Rogadinae (ER Insecta: AWH Hymenoptera: HERF} Braconidae) 
BET T 4p T RELIÉ UITL. BUM 95% ~ 100% CH RMR HAS BERUEEIZ8 DNA HY TES 10 FPA REAPS AL 5 PRON BE 
#PARAY 28S cDNA D2 Fr JEFE. (GenBank JF AY167645 - AY167659), FI BLAST 8 BAB AAAI, KAT GenBank 
13 1 FE2SR 28S rRNA D2 ARPS, RERIT TOT. FI 3 TOME GE S ARR) 3 ERIE GER 
distance-based neighbor-joining, NJ; E f £915; maximum parsimony, MP; AER DATE maximum likelihood; ML) 4rd J A Se ee PF 
BANATRARRAKA. GRR, ATR EBIASIRSIBU AT: 38 SUPE D] oA ee eS PAR. AARRE 
AKT BAER EA at GMER) RERRR AAA AAR CH. NJ. MP A ML oT E RER AH P Se 
Rogadini DÆ —4 HAR, ME- -THR BRI ERU RERNA. ARE IATER Top “A RETR Aleiodes + 
Hk Sa WE JR Arcaleiodes” W “HARRER Canalirogas + FEAT RIER Conspinaria + ABI Spinaria + AAIE Rogas”, — 7S 
ÆRE ARREA Aleiodes #15 Hki 8 JR Arcaleiodes BAX EAGT. FREI Aleiodes J& — TPA, FA AT FX 2 PMR 
Xo HRA EWA EAL A. SRRA Arcaleiodes Chen et He» 1991 J& — P REB. 4p sc IAA 
HRA Canalirogas + Ei REJA Conspinaria + i| BRE Spinaria + ARERR Rogas” BUB RTEDUG-ANIBEAT STR RHE: A 
Fa ARRP MFA RRA BR WER Spinaria: 42r T 2T ATI H UESCOX — Ex. MAREI Clinocentrini ÆA% 
A FEARBVAWARKRA PETES OR) Wie. RNMRARERM, PERE Yelicones UB MR ME [8 
Stiropius FAT IE A) RAEN EE B Y eliconini EP S EET. Stiropüni A 2 PZS, SEAN EDSHAR-R, BENEA 
TREEOIDURLES ZA SUAR. EE DE EESBR, CUR PE Ja SE DH. 

Aid: HESSE. Re; AREF; 28S rRNA: PTARRA 
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